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1 Introduction 


Interest for supersymmetric relics leans on the following properties: i) if R-parity is 
conserved, the Lightest Supersymmetric Particle (LSP) is stable; ii) if colourless and 
uncharged, the LSP is a nice realization of a relic Weakly Interacting Massive Particle 
(WIMP) [^. In the present paper we assume that supersymmetry exists in Nature and 
that properties (i) and (ii) hold. Thus the LSP will be our relic WIMP candidate; actually, 
after some preliminary considerations, our attention will be mainly devoted to the case 
when the LSP is provided by the neutralino. 

In Sect. 2 we discuss some generic connections between WIMP relic abundance and 
event rates for direct and indirect WIMP detection. Sect. 3 is devoted to a presentation 
of supersymmetric schemes. Numerical results and conclusions are hnally presented in 
Sect. 4. 

2 Relic abundance and detection rates 

When we investigate a WIMP as a relic in the Universe, the two most crucial questions 
are: a) does it contribute in a signihcant way to the matter content in the Universe , b) 
is it detectable either directly or through indirect signals? 

For a given WIMP, its cosmological relevance and its detectability are related, but in 
a way which is somewhat different from what one might naively believe. Actually, it turns 
out that, in the framework which is usually employed for the WIMP decoupling from the 
primordial plasma, WIMPs with highest detection rates might have a relatively modest 
relic abundance. For a detailed discussion of these features we refer to Ref. [^; here we 
simply sketch the main points. 

For dehniteness, our discussion is formulated in terms of relic neutralinos, though 
many features apply to a larger class of WIMP candidates, provided they share with the 
neutralino the property of having an elastic scalar cross-section off nuclei dominated by 
a coherent contribution in the range of relevance for current experiments. 

Under this circumstance, the detection rate R in WIMP direct measurements is pro¬ 
portional to the product Px'‘^scaiar^°^^ where is the neutralino local density and 
is the neutralino-nucleon cross-section. 

The values to be assigned to the neutralino local density = ^-pi have to be consistent 
with the values of the neutralino relic abundance When > {Qrnh‘^)min, where 

min is the minimum value of compatible with halo properties, we simply set 

Px ~ Pi When < {flmh‘^)min, the neutralino cannot be the unique 

cold dark matter particle, thus we assign to the neutralino a rescaled local density p^ = 

Pi X / {Vtmh'^)min (l-G., ^ / {Vtmh^)niin) 0 • 

Taking into account these rescaling properties of the local density, we hnd that R 
behaves as follows 
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where aann and v are the neutralino-neutralino annihilation cross-section and the relative 
velocity, respectively. < aann v > denotes the thermal average of the product {aann ■ v) 
integrated from the freeze-out temperature to the present-day one. In writing Eqs. (Ih 
H), we used the fact that in the standard scheme for the decoupling of a WIMP from the 
primordial plasma one has oc < aann v 

Now, the cross sections and aann, as functions of any generic coupling parame¬ 

ter C, behave similarly (he. they usually both decrease or increase in terms of variations of 
this parameter), because of crossing symmetry. Thus, for instance, ^ann are 

both increasing functions of tan (3 (see later on for dehnition), when the relevant processes 
are mediated by Higgs bosons. Usually, increases somewhat faster than aann, or 

approximately at the same rate. Thus, due to the properties displayed in Eqs. the 

typical behaviour of the rate R may be summarized as follows: i) for small values of (, 
both aud aann are small (and then is large) and R grows proportionally to 

^scaiar °’^^, h) as the strength of the coupling increases, in the region beyond the value (r (at 
which = {Vtmh?)min), Eq. (@) applies: the rate R still increases (though less rapidly 
than before rescaling), or remains approximately flat. This behaviour obtains up to values 
of the relic densities which fall short of the cosmological interesting range by a couple of 
decades. The features just discussed for the WIMP direct detection also apply to the 
indirect experiments consisting in detecting a neutrino flux arriving from the Sun or from 
the center of the Earth, due to possible WIMP pair annihilations in these two celestial 
bodies. Indeed, the size of these neutrino fluxes depends on the product . 

On the basis of the previous arguments, we can conclude that, in the cases of di¬ 
rect detection and indirect detection through pair annihilation in celestial bodies, the 
detectability of relic neutralinos is usually favoured for neutralinos of small i. e. 

for neutralinos which comprise only a subdominant dark matter component. 

Therefore, when one explores the possibility of detecting a specihc kind of WIMP, 
it would be self-defeating to restrict the analysis to the case when the putative WIMP 
contributes dominantly to the matter content in the Universe. Because of this fact, at 
variance with many analyses by other authors, our scanning of the susy parameter space 
will also include conhgurations entailing small neutralino relic abundances. Obviously, 
among the different situations which one may hnd in the course of the investigation, 
the case when the WIMP candidate is both detectable and of cosmological interest will 
represent the most rewarding occurrence. In the last section, we show that, indeed, a 
number of neutralino conhgurations of cosmological interest, though disfavoured by the 
previous arguments, may reach the level of detectability by current experiments. 
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In Ref. 0 it is shown that for processes depending on pair-annihilation in the halo the 
maximal rates occur for values of the relic abundance around the value Sub¬ 

dominant neutralinos are disfavoured for detectability by this type of signals as compared 
to neutralinos with a relic abundance around the value ~ 


3 Which model for supersymmetry? 

Supersymmetry, though strongly motivated by theoretical arguments, is not yet sustained 
by experimental evidence. Only some possible hints are available: unihcation of the 
Standard Model (SM) coupling constants at a Grand Unihcation (GUT) scale [Q, Higgs 
events at LEP2 i, I, 0, i i, ht of precision electroweak data improved by susy effects 

Eg. 

However, in the following we assume that supersymmetry exists in Nature and that 
R-parity is conserved (thus the LSP is stable). The nature of the LSP depends on the susy- 
breaking mechanism and on the specihc regions of the susy parameter space. We consider 
here gravity-mediated schemes, and domains of the parameter space where the LSP is 
the neutralino. Extensive calculations on relic neutralino phenomenology in gravity- 
mediated models have been performed (among the most recent references see, for instance, 

MM)- 

Here we refer to an analysis which we performed in the Minimal Supersymmetric ex¬ 
tension of the Standard Model (MSSM) in a variety of different schemes, from those based 
on universal or non-universal supergravity, with susy parameters dehned at the grand uni¬ 
hcation scale (GUT), to an effective supersymmetric model dehned at the Electro-Weak 
(EW) scale. 

3.1 Universal and non—universal SUGRA 

The essential elements of the MSSM are described by a Yang-Mills Lagrangian, the 
superpotential, which contains all the Yukawa interactions between the standard and 
supersymmetric helds, and by the soft-breaking Lagrangian, which models the break¬ 
ing of supersymmetry. The Yukawa interactions are described by the parameters h, 
which are related to the masses of the standard fermions by the usual expressions, e.g., 
rrit = h^V 2 ,rnb = h^Vi, where Vi are the nen’s of the two Higgs helds. Hi and H 2 . Im¬ 
plementation of this model within a supergravity scheme leads naturally to a set of uni¬ 
hcation assumptions at a Grand Unihcation (GUT) scale, Mcut'- i) Unihcation of the 
gaugino masses: Mi^Maur) = ''^ 1 / 2 , h) Universality of the scalar masses with a common 
mass denoted by mo: mi^McuT) = hi) Universality of the trilinear scalar couplings: 
^\Mgut) = A'^{Mgut) = A'^{Mgut) = Homo. 

This scheme is denoted here as universal SUGRA (or simply SUGRA). The relevant 
parameters of the model at the electro-weak (EW) scale are obtained from their corre¬ 
sponding values at the Mgut scale by running these down according to the renormalization 


11, 12, 13, 14, 15 
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group equations (RGE). By requiring that the electroweak symmetry breaking is induced 
radiatively by the soft supersymmetry breaking, one hnally reduces the model parameters 


to hve: mi/ 2 , mo, Rq, tan/9(= V 2 /V 1 ) and sign fi. 

We remark that in this very strict scheme the phenomenology of relic neutralinos is 
very sensitive on the way in which various constraints (for instance those on the bot¬ 
tom quark mass, rrib, on the top quark mass, rrit, and on the strong coupling cts) are 
implemented. 

Models with unihcation conditions at the GUT scale represent an appealing scenario; 
however, some of the assumptions listed above, particularly ii) and hi), are not very solid, 
because, as was already emphasized some time ago |T8|, universality might occur at a 
scale higher than Mqut ~ 10^® GeV, e.g.^ at the Planck scale. 

An empirical way of taking into account the uncertainty in the unification scale consists 
in allowing deviations in the unification conditions at Mcut- For instance, deviations from 
universality in the scalar masses at Mcuti which split from may be parametrized 
as M^j.^Mcut) = '^ 0(1 + This is the case of non-universal SUGRA (nuSUGRA) that 
we considered in Refs. [^, 0. Further extensions of deviations from universality in 

SUGRA models which include squark and/or gaugino masses are discussed, for instance, 
in 
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More recently, the possibility that the initial scale for the RGE running, M/, might 
be smaller than Mgut ~ 10^® has been raised, on the basis of a number of string models 
(see, for instance, 0,13, 01 and references quoted therein). As is stressed in Ref. P3| , 
Ml might be anywhere between the EW scale and the Planck scale, with signihcant 
consequences for the size of the neutralino-nucleon cross section. 


3.2 Effective MSSM 

The large uncertainties involved in the choice of the scale Mj make the SUGRA schemes 
somewhat problematic: the originally appealing feature of a universal SUGRA with few 
parameters fails, because of the need to take into consideration the variability of Mi 
or, alternatively, to add new parameters which quantify the various deviation effects 
from universality at the GUT scale. Thus, it appears more convenient to work with a 
phenomenological susy model whose parameters are defined directly at the electroweak 
scale. We denote here this effective scheme of MSSM by effMSSM. This provides, at 
the EW scale, a model, defined in terms of a minimum number of parameters: only 
those necessary to shape the essentials of the theoretical structure of an MSSM, and of 
its particle content. Once all experimental and theoretical constraints are implemented 
in this effMSSM model, one may investigate its compatibility with specihc theoretical 
schemes at the desired M/. 

In the effMSSM scheme we consider here, we impose a set of assumptions at the 
electroweak scale: a) all trilinear parameters are set to zero except those of the third 
family, which are unified to a common value A; b) all squark soft-mass parameters are 
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taken degenerate: rriq. = rriq] c) all slepton soft-mass parameters are taken degenerate: 

= mp d) the U{1) and SU{2) gaugino masses, Mi and M 2 , are assumed to be linked 
by the usual relation Mi = (5/3) tan^ 9wM2 (this is the only GUT-induced relation we are 
using, since gaugino mass unihcation appears to be better motivated than scalar masses 
universality). As a consequence, the supersymmetric parameter space consists of seven 
independent parameters. We choose them to be: M 2 , p., tan/?, rUg, mp A and vary 
these parameters in the following ranges: 50 GeV < M 2 < 1 TeV, 50 GeV < |p| < 
1 TeV, 80 GeV < tua < 1 TeV, 100 GeV < rriq, mj < 1 TeV, —3 < A < +3, 1 < 
tan/3 < 55 {rriA is the mass of the GP-odd neutral Higgs boson). 

The effMSSM scheme proves very manageable for the susy phenomenology at the 
EW scale; as such, it has been frequently used in the literature in connection with relic 
neutralinos (often with the further assumption of slepton/squark mass degeneracy: rriq = 
mj). Notice that here we are not assuming slepton/squark mass degeneracy. 

We recall that even much larger extensions of the supersymmetric models could be 
envisaged: for instance, non-unihcation of the gaugino masses ^ 

GP-violating phases p3[ . 


and schemes with 


Here we only report results in the effective scheme at EW scale, except for a few 
comments on universal and non-universal SUGRA results in Sect. 4. For further details 
we refer to Refs. am. 

The neutralino is dehned as the lowest-mass linear superposition of photino ( 7 ), zino 
{Z) and the two higgsino states H^)'. x = ®i7 + 0 . 2 Z + a^Hl + a^H^. Hereafter, the 
nature of the neutralino is classihed in terms of a parameter P, dehned as P = a^+a^- The 
neutralino is called a gaugino when P > 0.9, a higgsino when P < 0.1, mixed otherwise. 


3.3 Constraints on supersymmetric parameters 

In our exploration of the susy parameter space, we have implemented the following ex¬ 
perimental constraints: accelerators data on supersymmetric and Higgs boson searches 
(GERN collider LEP2 [Q and Gollider Detector GDF at Fermilab |^); measure¬ 
ments of the b ^ s + x decay . 

The new measurement of the muon anomalous magnetic moment was recently 

considered by a number of authors as a sign of supersymmetry. However, a deviation 
of the experimental value from the SM evaluation never exceeded a mere 1.8 a, or 


even less ||^, due the uncertainties affecting the evaluation of the hadronic contributions. 
Now, in the light of the recent clarihcation about the sign of the hadronic light-by-light 
contribution , the putative deviation effect on is evaporating away. Thus, in our 

analysis we do not assume that the new determination of is a sign of supersymmerty; 
we only use as a constraint: the supersymmetric contribution to muon anomalous 
magnetic moment is set to be constrained in the range: —200 < 0 ®“*^ ■ 10^^ < 640. 
To establish this interval, we have combined the results of Refs, 
hadronic vacuum polarization contributions and Refs. 


31, 32 


for the 


29, 301 for the hadronic light-by- 
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light contributions. 

As we mentioned in the Introduction, we have not restricted our exploration of the 
supersymmetric parameter space by requiring the neutralino relic abundance VL^h'^ to sit in 
any a priori chosen cosmological range. This will enable us to analyse our results in terms 
of relic neutralinos comprising either a dominant or a sub-dominant relic population. 


4 Results and conclusions 


In what follows we present some of our results in the framework of the effMSSM scheme. 


The evaluation of follows the procedure given in 


The cross-section has been calculated with the formulae reported in Refs. 


Important entries in crscaiar°”^ ^.re the quantities < N\qq\N >, where the quark scalar 
densities qq are averaged over the nucleonic state. The values of these quantities, derived 
from the pion-nucleon sigma term a-jrN and other hadronic quantities, are affected by 
large uncertainties |^. The quantity < A^|ss|A^ >, which is the most important term 
among the < N\qq\N >’s unless tan/? is very small |^, is affected by an uncertainty 
factor larger than 3. This conclusion is reinforced by the most recent determinations of 

H- 

3 employ the following set of values for the quantities 
^]): mi < U > = 23 MeV, mg < ss > = 


o'ttn HI; discussed in Ref. 

The results presented in Figs. 1 
rUg < N\qq\N > (denoted as set 1 in Ref. 

215 MeV, < hh > = 50 MeV. For the reasons discussed above, all values concerning 
in Figs. 1-3 are subject to an increase of ~ 4, when the current uncertainties in 


cr. 


(nucleon) 
scalar 


the mq < N\qq\N >’s are taken into account. 

Now let us turn to the presentation of our main results. In Fig.l we give the scatter 
plot for versus The two vertical lines denote the favorite range for VLmh?'- 

0.05 < < 0.3. The two horizontal lines bracket the range of sensitivity in current 


WIMP direct experiments [^, which, taking into account astrophysical uncertainties 
turns out to be 4 ■ 10“^° nbarn < < 2 ■ 10“® nbarn, for WIMP masses in 

the interval 40 GeV < myy < 200 GeV. Fig.l shows that the present experimental sen¬ 
sitivity in WIMP direct searches allows the exploration of supersymmetric conhgurations 
compatible with current accelerator bounds; a number of conhgurations stay inside the 
region of cosmological interest. 

Once a measurement of the quantity is performed, values for the local 

density versus the relic abundance may be deduced by proceeding in the following 
way 1) p^ is evaluated as [Px'C^scaiar Uxpt / cr/,aiar 4 where [p^-f^scaiar Uxpt denotes 
the experimental value, and is calculated as indicated above; 2) to each value 

of p^ one associates the corresponding calculated value of The scatter plot in 

Fig.2 is derived for the representative value [p^/(0.3 GeV cm“^) ■o'scaiar°°^]ea:pt = 1 ■ 10“® 
nbarn within the annual-modulation region of Ref. [^, and by taking in the range 
40 GeV < mw < 200 GeV. 


6 


































effMSSM 


C 


c 

Zu 

Z£ 




10 -» 

10 -’ 

10 -“ 

10 -" 

10-lD 

10 -" 

10 -'® 

10-ia 

10 -'^ 
10 -“ 


A. Boitinv. N. Forntngo. S. Scop«l (SOOI) 


10 -' 


n 


Figure 1: Scatter plot of cTgcaiar”"^ versus is taken in the the interval 40 GeV < 

mw < 200 GeV. The two horizontal lines bracket the sensitivity in current WIMP 
direct experiments 4 • 10“^° nbarn < < 2 ■ 10“® nbarn, for ^ = 1. The two 

vertical lines denote the range 0.05 < < 0.3. Dots (crosses) denote gaugino (mixed) 

conhgurations. 


The plot of Fig.2 shows that the most interesting region, he. the one with 0.2 GeV 
cm”^ < Px — GeV cm”^ and 0.05 < < 0.3 (cross-hatched region in the hgure), 

is covered by susy conhgurations probed by the WIMP direct detection. Let us examine 
the various sectors of Fig.2. Gonhgurations above the upper horizontal line are incom¬ 
patible with the upper limit on the local density of dark matter in our Galaxy and must 
be disregarded. Gonhgurations above the upper slanted dot-dashed line and below the 
upper horizontal solid line would imply a stronger clustering of neutralinos in our halo as 
compared to their average distribution in the Universe. This situation may be considered 
unlikely, since in this case neutralinos could fulhll the experimental range for but they 
would contribute only a small fraction to the cosmological cold dark matter content. For 
conhgurations which fall inside the band delimited by the slanted dot-dashed lines and 
simply-hatched in the hgure, the neutralino would provide only a fraction of the cold 
dark matter at the level of local density and of the average relic abundance, a situation 
which would be possible, for instance, if the neutralino is not the unique cold dark matter 
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effMSSM 



Figure 2: Scatter plot of versus This plot is derived from the experimental 

value [px/(0.3 GeV cm“^) = 1 ■ 10“® nbarn and by taking in the the 

interval 40 GeV < mw < 200 GeV, according to the procedure outlined in the text. 
The two horizontal lines delimit the range 0.2 GeV cm“^ < Px — GeV cm“^; the 
two vertical ones delimit the range 0.05 < < 0.3. The band delimited by the two 

slanted dot-dashed lines and simply hatched is the region where rescaling of pi applies. 
Dots denote gauginos, circles denote higgsinos and crosses denote mixed conhgurations. 


particle component. To neutralinos belonging to these conhgurations one should assign a 
rescaled local density. 


We remind that the scatter plot in Fig.2 refers to a representative value of [py. 


■O'. 


(nucleon) 1 
scalar J 

inside the current experimental sensitivity region, thus the plot in Fig.2 shows that current 
experiments of WIMP direct detection are probing relic neutralinos which may reach 
values of cosmological interest, but also neutralinos whose local and cosmological densities 
may provide only a very small fraction of these densities. 

versus The solid 

when only the 


In Fig.3 we give the scatter plot for the quantity 


line denotes the frontier of the 3(T annual-modulation region of Ref. 141 


uncertainties in pi and in the dispersion velocity of a Maxwell-Boltzmann distribution, 
but not the ones in other astrophysical quantities, are taken into account. Effects due 
to a possible bulk rotation of the dark halo or to an asymmetry in the WIMP velocity 
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Figure 3: Scatter plot of versus Crosses (dots) denote configurations with 

> 0.05 < 0.05). The solid contour denotes the 3 ct annual-modulation region 

of Ref.|^^ (with the specihcations given in the text). 


Our results in 


distribution would move this boundary towards higher values of | 

Fig.3 show that the susy scatter plot reaches up the annual-modulation region of Ref.E 


In our hgures only results referring to the effMSSM scheme are reported. For compar¬ 
isons among results in various schemes: universal SUGRA, non-universal SUGRA and 
effMSSM, we refer to Refs. [§, ]R|. As for the universal SUGRA, we only wish to remark 


that this very constrained model, combined with the present rather stringent experimen¬ 
tal bounds from LEP2, typically entails a sizeable suppression of the neutralino-nucleon 
cross-section. Whether or not this suppression may prevent the calculated o'scaiar°’^^ from 
reaching the region of present experimental sensitivity does depend on how the various 
constraints (typically the bounds on Higgs masses, on rrit and rrih) are implemented in 
the evaluations. By way of example, it is worth mentioning that the explicit bounds on 
the quantity sin^(Q; — /?) (a being the Higgs mixing angle in the neutral CP-even Higgs 
sector) as a function of ruh should be taken into account, rather than using a flat lower 
bound of 115 GeV for mh P, ||. Neglecting these features in a SUGRA calculation may 
lead to biased conclusions. 

We now summarize the main points of this paper: 
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• Most recent theoretical developments suggest supersymmetric schemes which no¬ 
tably differ from a strict model such as the universal SUGRA and point to the fact 
that this constrained scheme should be relaxed in many instances. Here, we have 
employed an effective MSSM scheme at the electroweak scale, which is particularly 
convenient to treat the relic neutralino phenomenology. 

• We have shown that current direct experiments for WIMPs, when interpreted in 
terms of relic neutralinos, are indeed probing regions of supersymmetric parameter 
space compatible with all present experimental bounds. 

• We have proved that part of the configurations probed by current WIMP experi¬ 
ments entail relic neutralinos of cosmological interest, and, a fortiori also neutralinos 
which might comprise only a fraction of the required amount of dark matter in the 
Universe. 
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